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Introduction 

This  paper  deals  with  the  preparation  of  potential  energy 
curves  for  diatomic  molecules.  The  starting  point  for  the 
calculations  is  the  spectroscopic  data  available  which  can 
include  observations  of  transitions  between  rotational, 
vibrational,  and  electronic  states.  Once  this  data  is  gathered, 
how  is  it  processed  so  that  the  entire  potential  energy  curve 
can  be  calculated. 
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Modeling  Spectroscopic  Data  for  Diatomic  Molecules 


Spectroscopic  data  resulting  from  the  emission  of  energy  by  a 
diatomic  molecule  during  a  transition  from  an  excited  electronic 
state  to  a  lower  state  may  be  represented  as  follows: 


v(v* ,J* ,v", J")  =  T'(v'fJ')  -  TM(v",J«) 


where  v  is  the  observed  line  frequency  in  wave  numbers,  the  term 
values  for  the  upper  and  lower  electronic  states,  T'  and  T" 
respectively,  are  represented  by  the  Dunham  type  expression^ 1 ^ : 

T(v, J)  s  A^Cv  +  1/2)1  Jj  (J  +  1)j  (2) 

Although  Eq  (2)  has  the  same  form  as  that  associated  with  the 
Dunham  coefficients,  Y.  i.e., 


>,J)  s  1I  Yij(v  +  1/2)1  Jj  (J  +  1)J 


researchers  M.  M.  Hessel  and  C.  R.  Vidal  make  the  distinction 
that  their  A^j  values  are  not  necessarily  identical  to  Dunham's 

Yij’s- 

(2) 

Dunham  arrived  at  the  expression  in  Eq  (2)  by  using  the 
Wentzel-Kramer-Brillouin  method  (WKB)  to  solve  the  Schroedinger 
Equation  for  the  rotating  diatomic  vibrator  for  the  energy  levels 
within  one  electronic  state: 


d2\p  8ir2yr2 

— »  +  - E  -  V(?)  - 


— *  + 

dr 


h  j ( J+1 ) 


8ir*r“  V(1+€)' 


ip  S  0 


•  W  "j  v  v  A.  'j-  ->> 


r  ••  - ",  ■  ■  v 


where: 


5  =  (r  -  re)/re 

r@  =  the  equilibrium  internuclear  separation 

u  a  the  reduced  mass 

V  =  the  potential  of  the  function,  with  a  minimum  at  rg 

Vr  =  the  last  term  which  is  due  to  the  centrifugal  force  of 

rotation . 

He  expressed  the  potential  energy  by  expanding  V  about  5  =  0. 

V  =  hcaQ52(1  +  a1  5  +  a2  52  +  a3  53  +  ...)  (5) 

2 

where  ag  =  u>e  /4Bg;  we  is  the  classical  frequency  of  small 

2  o 

oscillations  expressed  in  cm-1  and  =  h/(8Tt  \ir~  c).  Finally, 

6  6 

he  obtained  expressions  for  the  Y^'s  in  terms  of  the  a^'s  in  Eq 
(5). 

Hessel  and  Vidal  distinguish  their  constants,  A^,  from  the 

Dunham  constants,  Y^,  by  pointing  out  .that  the  latter  are 

theoretically  derived,  while  the  A^j's^aVe  the  results  of  a 

least-squares  fit.  Inherent  in  the  least-squares  fitting 

procedure  is  the  fact  that  the  A^  coefficients  absorb  the 

effects  of  inaccuracies  in  the  data,  and  are  somewhat  dependent 

upon  the  values  of  their  neighboring  constants  and  upon  any 

(•>) 

missing  constants  .  For  these  reasons,  the  A^’s  are  only 
estimates  of  the  theoretical  Y^  derived  by  Dunham. 

Some  further  discussion  of  the  Ajj  and  constants  is 
appropriate.  First,  a  listing  of  the  correspondence  of  Dunham’s 


constants  and  the  classical  spectroscopic  constants  is  desirable 
They  are  given  in  Table  1 . 


TABLE  1 


This  is  the  expression  normally  reported  in  the  literature  based 

(2) 

upon  the  work  of  Dunham  ,  although  the  expression  for  Yqq  is 

(7 ) 

not  explicitly  stated  in  Dunham's  work.  Sandeman  '  and 
( 8 ) 

Jarmain  expanded  on  Dunham's  work.  Using  this  material,  the 
expression  given  in  Eq  (6)  can  be  verified.  In  addition  to  this 
derivation,  the  expression  for  YQ0  was  found  in  the  following 
works: 10* 1 1 ' ^ .  The  expression  for  YqQ  was  found  to  be 
somewhat  different  in  two  publications.  First,  Herzberg's  YnA  i 


larger  than  in  Eq  (6)  by  +3/4  Ygg^1^.  Second»  in  McKeever's 
work^11^  due  to  a  typographical  error,  "  +  1/4  Y20"  was  omitted. 

The  works  of  Sandeman  and  Jarmain  compliment  Dunham's  work, 
making  easier  the  understanding  of  Dunham's  work. 

Returning  to  Eq  (1),  it  should  be  noted  that  when  a  least- 
squares  fit  is  performed  using  Eq  (2),  the  bottom  of  the 
potential  energy  curve  of  the  lower  state  is  typically  assigned  a 

ft 

value  of  ,,0W .  This  is  done  by  omitting  the  AQ0  term  from  the 

t 

fit.  When  Ann  is  found  from  the  least-squares  fit,  it  is  not 


equivalent  to  Dunham  's  Y00  for  the  upper  state.  For  this  reason, 

*  * 

Vidal  assigns  an  asterisk  to  Ann  and  calls  it  Ann.  This  term  is 


made  up  of  three  components  as  shown 


(14). 


Aoo  s  A00  +  Te  “  Aoo 


I  II  I  II 

Ag0  and  Aqq  are  estimated  from  Eq  (6).  This  Aqq  and  AQq  are  then 

f  t! 

the  best  estimates  of  Dunham's  Yqq  and  Ygg.Eq  (7)  can  be  solved 
for  Te',  the  electronic  term  energy  of  the  upper  state: 


»  I  w 

Te  =  A00  ”  A00  +  A00 


In  applying  Dunham's  expression,  one  should  be  aware  of  its 

(10) 

limitations.  Expressions  such  as  the  Morse  potential  : 


U ( r-r)  =  D  (1  -  exp(-B(r-r  ) )' 
e  e  e 


5 


are  constructed  to  guarantee  that  as  r  +  rg)  U  -*•  0  and  as  r  -*•  00 , 

U  approaches  the  dissociation  energy,  Dg.  Eq  (5)  on  the  other 
hand  does  not  necessarily  satisfy  these  two  criteria.  Sandeman 
showed  that  if  Morse's  equation  is  expanded  about  (r-r  )  =  0,  it 
takes  on  a  form  similar  to  Dunham's  Eq  (5),  but  is  expressed  in 
terms  of  two  constants  and  additional  numerical  coefficients  as 
follows : 

U  =  aQ(1  -  a+.583a2£-0.250a3£+0.086a4£-. . .)  (10) 

A  comparison  of  Eqs  (5)  and  (10)  shows  that  Dunham's  expression 

is  a  more  general  form  of  Morse's  equation.  Sandeman  also 

applied  a  similar  expansion  to  an  equation  for  potential  energy 

( 1  c  ) 

developed  by  Kratzner  J  and  obtained  a  similar  correspondence 

with  Dunham's  expression.  Because  of  the  increased  number  of 

variables  and,  in  turn,  the  increased  flexibility,  Dunham's 

expression  has  the  potential  for  being  more  accurate.  But, 

because  of  this  freedom,  obtaining  constants  which  satisfy  the 

convergence  criterion  is  more  difficult  and  is  not  guaranteed  for 

large  values  of  r  or  v. 

( 8 ) 

Jarmain  made  a  term-by-term  comparison  of  similar 
expressions  developed  from  Dunham's  work  and  by  the  RKR  method 
and  3howed  that,  upon  neglecting  Dunham's  small  corrections,  the 
potentials  produced  by  the  two  are  mathematically  identical. 

Based  upon  this,  the  use  of  Dunham's  expression  (Eq  (3))  to 
represent  spectroscopic  data  for  input  to  an  RKR  program  for 
generating  potential  curves  is  justified  and  appropriate. 


In  applying  Dunham’s  expression,  one  should  remember  that  his 
equations  are  developed  for  small  oscillations  about  r  .  At 
energy  levels  approaching  the  dissociation  limit,  Eq  (2)  may  no 
longer  be  appropriate.  Eqs  (1)  and  (2)  are  for  simple  cases.  In 
cases  where  A-type  doubling  occurs,  or  where  isotopic  effects  are 
being  considered,  more  complex  models  must  be  substituted  for  Eqs 
(1)  and  (2).  Works  cited  in  the  bibliography  by  Vidal  and  his 
coworkers  present  variations  of  Eq  (2). 

Spectroscopic  Constants  by  Merging  Least-Squares  Fit  Data 
(a)  Development  of  the  Method 

The  large  quantity  of  data  obtained  in  a  spectroscopic 
analysis  requires  that  the  data  be  reduced  to  a  more  manageable 
form.  For  example,  Eqs  (1)  and  (2)  in  the  previous  section  are 
used  to  calculate  reportable  constants.  To  make  these  models 
usable,  the  A^  coefficients  must  be  determined.  The  most  widely 
used  method  of  generating  these  coefficients  is  to  perform  a 
linear  least-squares  fit  to  the  experimental  data.  For 
spectroscopic  data,  a  model  is  created,  typically  involving  a 
power  series  expansion  based  upon  the  rotational  and  vibrational 
quantum  numbers  as  expressed  in  Eqs  (1)  and  (2).  If  only  band 
head  data  is  being  analyzed,  the  model  takes  the  following  form: 


v  =  ii:0AiO  (v’  +  1/2)1  "  iE0AiO  (v"  +  1/2)1  (11) 

If  several  different  transitions  are  involved,  different 
values  of  the  A”  ground  state  constants  may  be  obtained.  In 


addition,  the  accuracy  of  the  A”'s  may  be  different.  It  is 
desirable  to  merge  these  constants  and  obtain  a  best  estimate  of 
the  ground  state  constants.  This  can  be  done  by  accomplishing  a 
weighted  least-squares  fit  of  the  data.  Such  a  merging  can 
produce  more  accurate  estimates  of  the  constants  and  reduce  the 
error  limits  associated  with  them.  The  merged  fit  may  also 
improve  the  upper  state  constants. 

Merging  may  also  be  desirable  for  the  case  when  two  groups  of 
data  are  available  with  significantly  different  standard  errors. 
For  example,  if  infrared  data  for  vibration-rotation  band 
transitions  and  microwave  data  for  transitions  between  adjacent 
rotational  levels  are  available,  the  greater  accuracy  of  the 
microwave  data  can  improve  the  accuracy  of  the  other  constants  if 
a  weighted  merged  fit  is  performed. 

There  are  other  methods  of  data  fitting  available.  These 

methods  include  least-squares  deviation,  maximum  likelihood,  and 

minimum  chi-squared.  All  have  different  fit  criteria^ .  Also, 

(17  181 

nonlinear  least-squares  fits  are  possible  ' '  .  The  present 

report  deals  only  with  the  least-squares  fit  techniques. 

The  least-squares  method  or  ’’regression”  method  minimizes  the 
sum  of  the  squared  deviations  between  observed  values  and  values 
calculated  using  the  constants  obtained  from  the  fit.  The  least- 
squares  method  provides  the  minimum  variance  linear  unbiased 
(MVLU)  estimates  of  the  constants;  that  is,  the  least-squares 
method  introduces  no  bias.  The  fits  produced  are  linear 
functions  of  the  data.  When  the  fit  is  used  to  reproduce  data, 


the  generated  data  exhibits  the  smallest  possible  variance  from 
the  original  data  that  can  be  achieved  with  the  model  and  data 
used(,6>. 

In  performing  a  least-squares  fit,  the  following  assumptions 
are  made: 

1.  The  model  (equation)  chosen  to  represent  the  data  is  a 
perfect  description  of  the  physical  event. 

2.  The  model  is  linear  in  the  constants  to  be  estimated. 

3.  The  mean  error  of  the  experimental  data  is  zero.  Any 
systematic  error  present  must  be  small  compared  to  the  variances 
and  random  errors  of  the  data. 

4.  The  variance-covariance  elements  must  be  finite  and  their 
relative  values  known  if  different  groups  of  data  are  to  be 
merged. 

In  reducing  spectroscopic  data,  three  methods  are  available: 

1.  The  reduction  of  each  band,  i.e.,  each  a  group  of 
transitions  from  one  electronic  state  to  another,  separately. 

2.  The  reduction  of  a  number  of  bands  simultaneously. 

3.  The  reduction  of  the  bands  separately  and  then  merging 
the  resulting  data. 

The  first  and  third  techniques  are  developed  in  this  paper. 
The  method  of  approach  is  now  described.  Least-squares  fitting 
techniques  are  based  upon  the  following  matrix  equation: 


where  v  is  a  column  matrix  containing  the  experimentally  observed 
line  frequencies  expressed  in  wavenumbers,  X  is  a  matrix  made  up 
of  the  (v  +  1/2)  and  the  J(J  +1)  terms  as  dictated  by  the  model 
in  Eqs  (1),  (2),  or  (11)  as  chosen.  B  is  column  matrix  made  up 
of  the  Aij  constants  which  are  to  be  calculated.  The  e  column 
matrix  contains  the  unknown  errors  associated  with  each  observed 
experimental  data  point. 

Eq  (12)  is  then  solved  for  g  as: 

8  =  (XT  X)"1  XTv  (13) 

The  estimated  variance  of  the  fit  is  expressed  as  follows: 

a2  =  ( v  -  X8)T(  v  -  XB)/fB  (14) 

where  fm  is  the  degree  of  freedom  of  the  calculation.  The  degree 
of  freedom,  fm,  is  equal  to  the  number  of  experimental  values  v 
used  in  the  fit  minus  the  number  of  constants  to  be  obtained  in 
the  8  matrix.  The  variance-covariance  matrix  for  the  constants 
obtained  in  Eq  (13)  is  calculated  by  the  following  relation: 

0  =  a2(XT  X)"1  (15) 

The  diagonal  elements  of  the  9  matrix  are  the  variances  of  the 
jn3tants  and  the  off  diagonal  elements  are  their  covariances, 
li.  correlation  between  the  calculated  constants  is  expressed  by 
the  correlation  coefficients  which  are  obtained  using  the 
variances  and  covariances  in  the  following  manner: 


X-state  with  significantly  smaller  errors  than  those  obtained  in 
separate  fits  of  their  data. 


The  discussion  of  the  b  state  data  was  reserved  until  last, 
because  they  were  obtained  by  making  reassignments  of  many 
transitions  that  were  published  in  previous  works. 

The  b-state  constants  listed  in  Table  2,  Column  I  are  based 

upon  the  six  assignments  marked  by  asterisks  in  Table  6.  The 

assignments  were  made  by  Glessner^  D '  and  Snyderv  '  based  upon 

the  works  of  Kurylo,  et.  al.^°^.  The  variance  of  the  constants 

in  Table  3  is  remarkably  small  considering  that  only  six  data 

points  with  five  constants  were  used  yielding  a  degree  of  freedom 

of  one.  A  review  of  Kurylo's  tentative  assignments^0^  and 

(2Q) 

tentative  assignments  made  by  Oldenborg,  Dickson,  and  Zare 
show  a  much  greater  difference  between  the  observed  and 
calculated  values  than  would  be  expected  from  the  constants  of 
Table  2. 

Using  the  b  state  constants  of  Table  2,  new  assignments  were 
made  for  the  experimental  data  which  had  been  marked  as  belonging 
to  the  b-X  transition  by  Kurylo  and  Oldenborg.  It  was  possible 
to  make  new  assignments  for  29  of  34  observed  values  which  had 
been  reported  as  belonging  to  the  b  state.  Twenty-six  of  those 
were  used  in  a  least-squares  fit  to  generate  new  b  state 
constants.  DUNCON  was  used  to  perform  separate  least-squares 
fits  to  the  data  of  Snyder  and  Glessner,  and  to  the  data  of 
Kurylo  and  Oldenborg.  Then  DUNCON  was  used  to  merge  the  results 
of  these  fits. 


increased  error  limit  reflects  a  bias  in  one  or  both  of  the  data 
sets. 

Column  III  of  Table  3  reflects  the  most  conservative  error 
estimate  of  the  three  columns  in  Table  3-  It  is  also  based  upon 
the  largest  quantity  of  data. 

Presented  in  Table  5  art  the  constants  obtained  by  merging 

the  experimental  data  of  Glessner  and  Snyder.  Presented  for 

comparison  are  results  published  by  other  researchers. 

The  D-state  constants  were  obtained  in  a  separate  merged  fit 

from  that  used  to  obtain  the  other  constants  because  of  their 

( 28 ) 

large  error.  The  values  reported  by  Blooraenthal  produced  a 
better  fit  for  his  data;  hence,  his  constants  are  the  preferred 
constants . 

The  A-state  constants  reported  in  Table  5  provide  a 

reasonable  agreement  with  those  reported  by  Linton  and 
( 27 ) 

Broida  .  The  reported  constants  are  based  upon  the  merging  of 
data  reported  by  Glessner  and  Snyder. 

The  small  a-state  constants  reported  in  Table  5  rival  those 
reported  by  Linton  and  Broida  in  accuracy.  Table  3  shows  that 
the  merged  data  of  Glessner  and  Snyder  produces  constants  with 
error  limits  of  magnitude  comparable  to  that  reported  by  Linton 
and  Broida. 

The  ground  state  (X)  constants  reported  in  Table  5  agree  well 
with  those  reported  by  Linton  and  Broida.  The  merging  of 
Glessner's  and  Snyder'3  data  resulted  in  constants  for  the 


Merged 

onstants 


tate 


Te  (cm-1) 


We  (cm-1) 


WeXe  (cm*"1) 


Is 


his  Study 

D 

30115.8(46.1) 

542.2(54.9) 

-1.99  (13.1) 

Glessner  & 

B 

22282.80(3.67) 

489.50(2.63) 

0.149(0.363) 

nyder,  Refs 

A0  + 

19856.22(3.31) 

447.88(3.50) 

1.013(0.6007) 

25)  &  (26)) 

b2 

16335.44(2.89) 

428.43(1.21) 

-0.99552(0.0983) 

a 

16021.51 (1.69) 

484.239(0.808) 

2.949(0.109) 

X 

721.905(0.620) 

3.5900(0.0683) 

inton  and 

A0+ 

19862.6(1.5) 

444.3(0.8) 

0.54(0.12) 

roida,  Ref 

al 

16024.9(1.45) 

481.5(0.7) 

2.45(0.07) 

27) 

XO+ 

720.97(0. 36) 

3.536(0.025) 

loomenthal , 

D 

30197.0 

530.6 

1.05 

lef  (28) 

B 

22289.8 

496.3 

2.  33 

X 

722.3 

3.73 

ildenborg , 

lickensen  & 

are.  Ref  (29) 

b 

16 379+43 01 

^"This  value  is  based  upon  estimate  by  Oldenborg,  Dickensen  and  Zare 

hat  the  b  state  is  350+430  cm-1  above  the  a  state. 

2  ** 

Based  upon  six  data  points  marked  with  asterick  in  Table  7. 

-•Errors  reported  in  parenthesis  are  one  standard  deviation. 


are  the  unmerged  and  merged  estimates  for  the  D-X  transitions. 

The  D-X  data  is  merged  with  other  data  in  a  separate  fit  to 

prevent  its  large  variance  from  degrading  the  constants  for  the 

other  transitions.  Finally,  the  merged  results  of  Glessner  and 

Snyder's  data,  are  compared  in  Table  5  with  the  best  previous 

published  values  available ^27 ’ ^ ^ . 

Examination  of  Table  2  reveals  values  of  ui  for  the  X-state 

e 

ranging  from  715.7  to  731.5  cm”1.  DUNCON  permits  the  merging  of 
the  data  listed  in  Table  2  in  a  weighted-correlated  least-squares 
fashion. 

The  results  of  the  merging  are  presented  in  Table  3.  Error 

limits  for  w  of  the  X-state  are  reduced  from  a  maximum  standard 
e 

error  of  +11.3  cm”1  to  +0.620  cm-1  when  the  data  in  Columns  I  and 
II  of  Table  2  are  merged.  The  results  are  presented  in  Column  I, 
Table  3. 

Column  II,  Table  3  contains  the  smallest  error  estimate  for 
the  constants  of  the  a-,  A-,  and  X-electronic  states.  The 
variances  for  the  merged  fits  of  Columns  I  and  II  are  similar, 
1.0248  and  1.008,  respectively.  But  when  the  data  in  Columns  I 
and  II  are  merged  to  yield  Column  III,  Table  3,  the  estimated 
variance  of  the  fit  and,  in  many  cases,  the  standard  errors  of 
the  individual  constants  increase.  Examination  of  (Columns  I  and 
II)  the  a-  and  X-states  reveals  terra  electronic  energies,  Tg  with 
error  limits  of  comparable  magnitude  which  do  not  overlap. 

Because  of  this,  when  merged,  the  error  limits  reflect  the  equal 
weighting  of  the  almost  significant  differences  and,  hence,  the 


TABLE  4 


Spectroscopic  Constants  for  D-X  Transition  of  PbO 


(cm“^ ) 


Electronic 

State 

Constant 

Unmerged 

Merged 

C322.31 

[4.19] 

D 

Te 

29897.8(67.4) 

30115.8(46.1) 

We' 

490.0(30.8) 

542.2(54.9) 

We’Xe 

-16.53(7.72) 

-1.99(13.1) 

X 

We" 

606.1(33.3) 

720.989(0.401 ) 

We"Xe" 

-8.32(3.53) 

■  —  ■ — .  J 

3.5312(0.0296) 

Errors  Reported  in  Parenthesis  are  One  Standard  Error 


Source  of  Experimental  Data 

Elec¬ 

tronic 

State 

Constant 

I 

Merging  of 
Glessner  &  Snyder 
Data 

II 

Merging  of 
Linton  &  Broida 
Data 

III 

Merging  of 
Glessner,  Snyder, 
Linton,  &  Broida 
Data 

a 

Te 

16021.51(1.69) 

16024.06(1.16) 

16023.16(1.80) 

• 

We' 

484.239(0.808) 

481.800(0.439) 

481.488(0.769) 

We ' Xe ' 

2.949(0.109) 

2.4595(0.0434) 

2.4339(0.0818) 

b 

Te 

16335.44(2.89) 

16331.83(4.77) 

We' 

428.43(1.21) 

428.34(2.27) 

We ' Xe ' 

-0.99552(0.0983) 

-1.014(0.186) 

A 

Te 

19856.22(3.31) 

19864.11(2.14) 

19860.04(3.13) 

We’ 

447.88  (3.50) 

443.39(1.12) 

445.32(1.83) 

• 

We ' Xe ' 

1.013(0.6007) 

0.446(0.145) 

0.655(0.249) 

B 

Te 

22282.80(3.67) 

22280.52(7.08) 

We' 

489.50(2.63) 

489.79(5.25) 

We ' Xe" 

0.149 (0.363) 

0.225(0.720) 

X 

We" 

721.905(0.620) 

720.736(0.244) 

720.890(0.395) 

We"Xe" 

3.5900(0.0683) 

3.5201 (0.0171) 

3.5316(0.0292) 

Variance 
of  Fit 

1.0248 

1.008 

4.076 

^•Column  I  is  the  result  of  merging  with  DUNCON  the  separate  fits  in 

Columns  I  and  II  of  Table  2. 

2 

Column  II  is  the  result  of  merging  the  separate  fits  in  Column  III 
of  Table  2. 

3Column  III  is  the  result  of  merging  all  separate  fits  in  Table  ? 
with  DUNCON. 

Errors  reported  in  parenthesis  are  one  standard  error. 
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TABLE  2  -  (Continued) 


ransition 


Source  of  Experimental  Data 

Constant 

I 

II 

III 

Glessner 

Synder 

Linton  4  Broida 

We'Xe  ' 

2.04(2.74) 

0.479(0.8632) 

0.443(0.146) 

We" 

715.7(11.3) 

724.69(2.93) 

720.933(0.374) 

We"Xe" 

2.80(1.58) 

4.012(0.390) 

3.5334(0.0239) 

[52.02] 

[64.633 

Te 

22290.9(10.3) 

22286.22(9.06) 

We' 

495.70(4.34) 

483.00(8.13) 

We'Xe  ' 

2.22(1.39) 

i 

• 

o 

*~N 

• 

- » 

We" 

731 .52(6.99) 

721.11(3.85) 

We"Xe” 

5.06(1.11) 

3.532(0.477) 

^he  quantities  in  brackets  are  the  variances  of  the  least-squares 
fits  (Eq  (14)).  Errors  reported  in  parenthesis  are  one  standard 
deviation.  Column  I  contains  least-squares  fit  to  data  produced  by 
Glessner  (Ref  19).  Column  II  constants  are  based  upon  data  produced  by 
Snyder  (Ref  42).  Column  II  constants  are  the  results  of  performing  least- 
squares  to  data  published  by  Linton  and  Broida  (Ref  33). 


2The  constants  for  b-x  were  obtained  by  combining  the  data  observed 
by  Glessner  and  Snyder. 
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Constants  from  Separate 
for  Each  Electronic  T 


ares  F 

t  cm  1  ' 


Constant 


We’Xe 


We"Xe 


e 

e" 

We’Xe 


o" 


I 

Glessner 

[10. 

03] 1 

16021. 

54(2. 

81) 

482. 

93(1. 

26) 

2. 

814  (0 

.170) 

719. 

36(2. 

19) 

3. 

067(0 

.455) 

[2. 

035] 

16336. 

43(6. 

08)  2 

428. 

98(1. 

46) 

-0. 

944(0 

.114) 

722. 

00(1. 

92) 

3. 

.554(0 

.193) 

[41. 

70] 

19846 

.0(25. 

8) 

1  453. 

5(18. 

5) 

Source  of  Experimental  Data 


II 

Synder 


[9.875] 

16021.68(2.59) 
484.51 (1.09) 
2.970(0.146) 
722.785(0.961) 
3.712(0.105) 


[33.07] 

19858.52(6.88) 

446.30(4.92) 


III 

Linton  &  Broida 


[6.310] 

16023.12(1.38) 

481.831(0.438) 

2.4590(0.0432) 

720.166(0.469) 

3.4600(0.0462) 


[4.817] 

19864.79(2.38) 

443.31(1.12) 


5.  The  standard  error  of  the  constants  and  the  degrees  of 
freedom  involved  in  the  calculations.  This  data  will  enable  the 


reader  to  establish  his  own  confidence  limits. 

6.  The  estimated  variance  of  the  variance,  cr  ,  may  be 

helpful.  An  estimate  of  this  value  for  large  samples  of  normally 

(24) 

distributed  data  may  be  obtained  by  the  formula 

oa  =  a/[2fm]1/2  (23) 

Other  data  which  may  be  appropriate  include  the  variance- 

covariance  matrix  of  the  final  merged  constants  and  the  results 

of  the  testing  to  determine  the  normality  of  the  data  sample. 

Application  to  the  Spectroscopic  Analysis  of  Lead  Oxide  (PbO) 

The  experimental  data  of  the  AFIT  investigators, 

Glessner^2^ ,  and  Snyder^2^  was  reduced  using  DUNCON  to  yield 

spectroscopic  constants  for  the  X,  a,  b,  A,  B,  and  D  electronic 

states  of  PbO.  In  addition,  the  reported  data  of  Linton  and 
(27) 

Broida  was  used  to  obtain  constants  for  the  X,  a,  and  A 
states.  The  data  of  Glessner  and  Snyder  was  reduced,  error 
estimates  were  obtained,  and  the  data  was  merged  using  DUNCON. 
Lead-oxide  data  produced  by  Linton  and  Broida  was  also  reduced 
using  DUNCON  and  the  A-X  and  a-X  data  were  merged.  Also,  the 
results  of  all  three  investigations  were  merged  using  DUNCON. 

The  results  of  reducing  the  data  of  each  investigator  for  each 
transition  separately  is  presented  in  Table  2.  This  data  is  then 
merged  in  several  combinations.  The  constants  resulting  from 
these  mergings  are  presented  in  Table  3«  Presented  in  Table  4 


the  extension  of  the  potential  curve  to  its  dissociation 
limits(22). 

Before  merging  data  from  separate  least-squares  fits,  a  check 
should  be  made  for  systematic  errors.  A  simple  check  is  to 
verify  that  for  a  given  set  of  confidence  limits,  say  95%,  the 
error  limits  for  the  common  constants  obtained  in  two  independent 
fits  overlap  each  other.  If  they  do  not,  then  it  is  probable 
that  one  set  contains  a  systematic  error.  A  second  more 
sophisticated  test  involves  the  computation  of  the  confidence 
limits  for  the  differences  for  pairs  of  corresponding  constants 
checking  to  see  if  they  include  zero^1*^.  Calculation  of  the 
confidence  limits  of  differences  involves  student's  t-factor  and 

(  0-3  ) 

methods  are  given  by  Bennet  and  Franklin  J  ;  and  Dixon  and 
Massey^20\  A  more  extensive  discussion  of  the  covariance- 
variance  data  and  how  it  may  be  used  is  presented  by 
Albritton(l6»19). 

Finally,  Albritton  provides  a  recommendation  as  to  what  data 
should  be  provided  in  a  report  on  the  results  of  a  spectroscopic 
study.  He  lists  the  following  as  essential  data: 

1.  The  observed  line  numbers  (with  assignments). 

2.  If  room  permits,  the  variance-covariance  matrix  used  to 
merge  the  data. 

3.  The  model  used  to  represent  the  experimental  data. 

4.  The  estimated  molecular  constants. 
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With  the  confidence  limits  for  the  estimated  constants 

properly  established,  one  can  then  determine  if  the  constants  are 
(21 ) 

significant'  .  For  example,  if  a  constant  has  a  value  of  3  and 

an  assigned  error  limit  of  +5,  the  range  in  which  the  true  value 

of  the  constant  might  be  found  is  -2  to  +8.  This  range  includes 

0;  hence,  the  constant  is  not  significant.  If  the  assigned 

1  /2 

confidence  limits,  g^+te^  includes  zero,  then  the  constant  is 
insignificant  and  may  normally  be  discarded.  Two  exceptions  to 
this  rule  may  be  encountered. 

The  first  exception  concerns  the  case  where  a  non-zero 
correlation  between  two  molecular  constants  exists.  If  non-zero 
correlation  coefficients  are  involved,  standard  errors  must  be 
calculated  using  the  most  general  formulation  involving  both 
variances  and  covariances.  If  large  correlations  exist,  the 
rounding  of  a  constant  to  the  number  of  significant  digits 
dictated  by  the  standard  error  of  the  constant,  may  result  in 
loss  of  information,  i.e.,  rounding  may  introduce  unnecessary 
inaccuracies  into  the  calculated  data^1^. 

Albritton  suggests  that  one  digit  beyond  the  "one  standard 
error  digit"  should  suffice  for  most  fits.  A  simple  check  is  to 
use  the  rounded  constants  and  see  that  they  reproduce  to  some 
desired  accuracy  the  values  calculated  using  the  unrounded 
constants  obtained  from  the  fit. 

The  second  exception  occurs,  where  it  may  prove  necessary  to 
retain  a  constant  when  not  justified  by  the  standard  error,  in 


As  in  all  analyses ,  the  data  should  be  checked  to  insure  that 


the  errors  predicted  by  the  least-squares  fit  follow  a  normal 

distribution  pattern.  Albritton  provides  a  fairly  complete 

discussion  of  the  checks  commonly  used^®). 

Having  established  that  the  errors  of  the  data  are  normally 

distributed,  one  should  then  construct  confidence  limits  for  the 

constants  based  upon  the  degrees  of  freedom  and  the  confidence 

one  wishes  to  have  that  the  true  value  is  within  the  assigned 

limits.  The  placing  of  too  much  reliance  upon  error  limits 

1/2 

specified  by  one  standard  deviation,  1x0^  >  should  be  avoided. 

Even  if  a  large  sample  of  data  is  taken,  one  can  only  be  about 

1/2 

68%  certain  that  the  "true  value"  lies  within  a  range  of  +1x0 • 

If  a  95%  confidence  level  is  desired  and  the  degrees  of  freedom, 

1/2 

fm,  is  greater  than  30,  one  must  assign  limits  of  +2x0 ^  •  The 
multiplier  of  is  student's  t-factor.  It  is  named  after  its 

originator,  W.  S.  Gosset  (1876-1937).  It  is  a  function  of  the 
degrees  of  freedom  of  the  calculation  and  the  desired  confidence 
level,  "l-o".  Dixon  and  Massey  (20)  give  a  tabulation  of 
student's  t-factor  as  a  function  of  the  degrees  of  freedom  and 
the  degree  of  confidence.  The  confidence  limits  for  the  constant 
are  then  expressed  as  follows: 

1/2 

Bi  +  t^m'  1-°t)0  ii  (22) 


Having  established  the  limits,  one  can  assume  that  the 

probability  of  the  true  value  of  falling  within  the  limit  +t0 
is  "l-o"  where  a  may  vary  between  0  and  1. 


.>.*.*.*.*.  \«  *  4 


•/y, 
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will  hold  the  desired  merged  constants.  The  y  matrix  will  often 

contain  several  values  for  the  same  constant.  Each  constant  will 

M 

be  represented  only  one  time  in  the  g  matrix. 

To  perform  a  weighted,  correlated  least-squares  fit,  a  matrix 
composed  of  the  variance-covariance  matrices  from  a  separate 
least-squares  fit  is  required. 

This  matrix  is  given  the  symbol  "61" .  It  is  structured  as 
follows: 


01 


0 


0 


(18) 


where  0.j ,  02,  and  0^  are  the  variance-covariance  matrices  of  Eq 
(15).  The  formula  required  to  obtain  the  merged  constants  is  a3 
follows: 


0M  =  (XT0l”^X)”^  XT0I-1y 


(19) 


The  estimated  variance  of  the  new  merged  fit  is  given  by: 

a2  =  (y  -  X6M)T0I_1  (y  -  X0M)/fm  (20) 

where  f  is  the  degree  of  freedom.  The  variance-covariance 
m 

matrix  of  the  merged  constants  is  obtained  by  the  following 
calculation : 

0M  =  a2(xT0I_1X)"1  (21) 


12 


The  subscripts  refer  to  the  two  constants  for  which  the 
coefficient  expresses  the  correlation. 

As  mentioned  earlier,  the  performance  of  separate  fits  to 
each  band  will  yield  several  values  for  the  constants  associated 
with  the  lower  electronic  state.  The  variances  for  different 
bands  may  also  differ  significantly.  If  the  variances  of  the 
bands  are  the  same,  the  data  can  be  combined  in  one  large  data 
group  and  reduced  simultaneously.  If  estimates  of  the  errors  are 
available,  a  weighted  fit  based  upon  the  error  may  and  should  be 
performed. 

Albritton  describes  an  approach,  the  correlated  least-squares 
fit,  which  is  more  general  than  the  weighted  least-squares 
fit(16,19).  Albritton  showed  that  one  need  not  use  all  raw  data 
at  one  time  to  obtain  a  merged  fit.  He  showed  that  a  weighted 
simultaneous  multiband  fit  and  a  merged  band-by-band  fit  are 
equivalent.  The  merged  fit  is  presented  in  this  paper. 

The  bands  or  other  groups  of  data  are  first  reduced  in  a 
least-squares  manner  as  previously  discussed. 

Then  the  data  obtained  in  the  initial  least-squares  fit  may 
be  expressed  in  the  following  manner: 

y  =  Xp11  +  A  (17) 

where  y  is  a  column  matrix  made  up  of  the  constants  obtained  from 
the  separate  least-squares  fits  and p  M  is  a  column  matrix  which 


TABLE  6 

Assignments  of  the  Vibrational  Bands  for  the  b-X  Transition  of 
PbO  Based  on  Molecular  Constants  Obtained  by  the  MVLU  Technique 


Observed  (cm“^) 

Calculated  (crn-^-) 

Obs-Calc 

(cm) 

Assignment 

Reference 

21643 

21659.8 

-16.8 

14,1 

7 

21200 

21207.6 

-7.7 

13,1 

7 

19869.8 

19860.2 

9.6 

10,1* 

Present  Work 

18515 

18526.3 

-11.3 

7,1 

2 

17488.4 

17483.1 

5.3 

3,0 

Present  Work 

17212 

17206.1 

5.9 

4,1 

2 

17108 

17120.5 

-12.5 

7,3 

2 

17033 

17047.5 

-14.5 

2,0 

2 

16915 

16938.0 

-23.0 

5,2 

2 

16773 

16769.0 

4.0 

3,1 

2 

16614 

16613.5 

-9.5 

1.0 

2 

16488 

16499.5 

-11.5 

4,2 

2 

16434 

16428.7 

5.3 

7,4 

2 

16297 

16292.4 

4.6 

17,11 

2 

16215 

16214.1 

0.9 

14,9 

2 

16066.8 

16062.4 

4.4 

3,2* 

Present  Work 

15773 

15762.0 

11.0 

13,9 

2 

15623 

15627.0 

-4.0 

2,2 

2 

15363 

15363.2 

-0.2 

3,3 

2 

15281.9 

15285.4 

-3.5 

9,7* 

Present  Work 

14687.8 

14671.4 

15.6 

3,4 

Present  Work 

14180.1 

14178.5 

1.6 

8,8 

Present  Work 

13367.3 

13369.4 

-2.1 

0,4 

Present  Work 

13117.5 

13117.4 

0.1 

1,5* 

Present  Work 

12440.1 

12440.4 

-0.3 

1,6* 

Present  Work 

11968.9 

11978.2 

-9.3 

3,8* 

Present  Work 

The  26  assignments  used  to  obtain  the  constants  are  listed  in 
Table  6.  The  resulting  b  state  and  X  state  constants  are 
presented  in  Table  7.  For  the  26  values  used  in  the  new  fit,  a 
standard  error  of  +10.5  cm  was  obtained.  For  the  15  values  which 
the  previous  investigators  reported  ••observed-calculated"  values, 
their  standard  error  was  +47. 4  cm. 

As  shown  by  this  example,  by  use  of  such  a  merging  technique, 
a  correlation  with  previous  work  can  sometimes  be  performed 
leading  to  more  reliable  band  assignments. 

Rydberg-Klein-Rees  Potential  Energy  Curves 

The  Rydberg-Klein-Rees  (RKR)  model  for  the  potential  energy 
curve  is  developed  in  this  section.  Rydberg  started  with  the 
expression  for  the  total  energy,  E,  for  a  rotating-vibrating 
diatomic  molecule: 

E  =  p2/2y  +  V(r)  (24) 

p 

The  p  /2y  term  in  Eq  (24)  is  the  kinetic  energy  of  the  system,  p 
in  turn,  is  the  linear  momentum.  The  second  term,  V(r)  is 
defined  as  follows: 

V(r)  =  U(r)  +  </r2  =  U(r)  +  P2/2yr2  (25) 


where  U(r)  is  the  potential  energy  of  the  system,  P  is  the 
angular  momentum  of  the  system,  and  y  is  the  reduced  mass. 


TABLE  7 

Spectroscopic  Constants  for  b-X 
Transition  for  Lead  Oxide 


Te 

We 

WeXe 

16325.1(11.2) 

430.99(2.47) 

721.41(4.35) 

-0.757(0.165) 

3.700(0.441) 

Based  on  this  energy  expression,  Rydberg  developed  a 
procedure  for  graphically  determining  the  classical  turning 

(  31  ) 

points  for  a  diatomic  molecule  .  Klein  modified  Rydberg’s 
derivation  so  that  the  turning  points  could  be  calculated 

(■3  O) 

numerically  J.  Reese  obtained  quadratic  and  cubic  analytic 

( -3-3  \ 

solutions  approximating  Klein's  formulas  . 

Klein's  mathematical  development  is  given  here.  Terms  used 
in  the  development  are  shown  in  Fig  1  as  an  aid  in  following  the 
presentation.  The  classical  turning  points  are  designed  r'  and 
r".  Klein  defined  the  turning  points  in  terms  of  quantities  r1 
and  r~.  For  his  analysis,  r^  and  are  measured  from  rg ,  r^ 
being  positive  and  negative.  Then  r*  and  r"  are  defined  as 
follows: 


r '  =  r.  +  r 
1  e 


r "  =  r  +  r  ~ 
e  d 


Then 


dr'  =  dr. 


dr”  =  dr. 


dr  =  dr^  or  drg 


Vlvlvlv 


Klein  designated  the  width  of  the  potential  curve  as  2f  and  from 
Eqs  (26)  and  (27),  it  follows  that: 


2f  =  (r1  +  rg)  -  (r2  +  rg)  =  r.j-r2  (3D 

The  energy  in  Eq  (24)  may  be  expressed  as  a  function  of  ”1", 
the  action  integral  for  a  rotating  vibrator  and  wtcw.  These  terms 
are  defined  as  follows: 

I  =  §Pdr  (32) 

«  =  P2/(2u)  (33) 

The  following  relationships  which  involve  the  period,  t,  and 
the  average  position  of  a  harmonic  oscillator,  l/r ,  are  known: 

1/t  =  3E/3I  (34) 

l/r2  =  OE/3<)  (35) 

P  =  t(3E/3<)  (35 ' ) 

Using  Eqs  (32)  and  (34)  and  solving  for  the  differential  of  E  and 

adding  the  results,  the  following  is  obtained: 

5E  =  (  1/t)  <51  +  ( p/r )  5<  (36) 
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Inverting  Eq  (32),  the  following  is  obtained: 

t  =  3E/3I  (37) 

Setting  Eq  (36)  equal  to  zero  and  solving  for  P,  the  following  is 
obtained: 

6  =  -  pi/ 3  k  (38) 

From  the  equation  for  the  momentum,  p  =  yv  =  y(dr/dt),  one  can, 
by  integrating  over  one  cycle,  obtain  an  expression  for  the 
period  of  vibration  for  a  harmonic  oscillator: 

x  =  §y(dr/p)  (39) 

Solving  Eq  (24),  the  equation  for  energy,  for  p  and  substituting 

into  Eq  (39),  the  period  of  oscillation  may  be  expressed  as: 

x=  (y/2)  §  (E-V(r))1/2dr  (40) 

The  integral  in  Eq  (40)  can  be  broken  into  two  parts,  one  to 

be  integrated  from  rg  to  r^,  and  the  other  from  rg  to  rg  (as 
depicted  in  Fig  1).  This  operation  yields: 

r*  r* 

t  =  2(u)1/2  if  1  dr/(E-V(r))1/2  -  /  e  dr/(E-V(r) ) 1/2]  (41) 


The  factor  of  two  is  introduced  to  account  for  one  vibration 
cycle  from  rg  through  r^  and  r2  and  back  to  rg  again. 

Rearranging  Eq  (41)  and  recombining  the  integrals,  one  obtains: 

t  =  (2y)1/2  Z1  (drrdr2)/(E-V(r))1/2  (42) 

r2 

Now  to  change  the  variables  of  integration,  Klein  defined  the 
following  terms: 

y  =  V(r )  -  V(re)  (43) 

and 

x  =  E  -  V(re)  (44) 

where  x  is  constant  for  a  given  E.  Then  r  may  be  expressed  as  a 
function  of  y  as  follows: 


r  r  r1 (y) ,  r  >  re  (45) 

r  =  r2(y),  r  <  rQ  (46) 

and  when  at  rQ,  r^(0)  =  r2(0)  =  r  .  Subtracting  Eq  (43)  from 
Eq  (44): 


x  -  y  =  E  -  V(r ) 


(47) 
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and  if  the  "f"  of  Eq  (31)  is  function  of  y,  then  as  defined 

r.j(y)  -  r2(y)  =  2f(y)  (48) 

then  1 


dr.|/dy  -  dr2/dy  =  2df(y)dy 


or 


(49) 


dr1  -  dr2  =  2f'(y)dy  (50) 

Using  Eqs  (45),  (46),  (47),  and  (50)  to  change  the  variable 
of  integration  and  integration  limits  in  Eq  (42),  Klein  obtained: 

T  =  2  (2U)1/2  /  1  2  f '  (y)  dy/(x-y)1/2  (51) 

*2 


From  Eq  (43)  when  r  s  r  ,  then  V(r)  =  V(r  )  and  y  a  0.  Thus,  the 
integral  in  Eq  (51)  may  be  evaluated  from  0  (defined  as  the 
bottom  of  the  potential  well)  to  some  value  X. 

t  s  2  (2u)1/2  /*  f ' (y)  dy/(x-y)1/2  (52) 

0 

This  is  the  first  of  two  integral  equations  central  to  Klein's 


derivation . 


Combining  Eqs  (28),  (35),  (35'),  and  (40)  yields: 

p  =  T  (1/r2)  =  (y/2)1/2  §  (1/r2)  [dr/ (E-V(r) ) 1 /2]  (53) 

By  definition: 

dr/r2  =  -(1/r)  (54) 

Then  Eq  (53)  becomes: 

p  =  -(y/2)5  §  d(1/r)/(x-y)1/2  (55) 


Again,  the  integal  in  Eq  (55)  may  be  broken  into  two  parts.  Then 

for  r  >  r  : 

e 

r"  s  r2  +  re  (56) 

1/r"  =  1/(r2  +  re)  (57) 

and  if  rg  is  set  equal  to  0  then: 

1/r"  =  1/r2  (58) 

Similarily  for  r  <  0: 

1/r'  =  1/r,  (59) 
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Then  for  small  oscillations  about  r  ,  Eq  (53)  may  be  written  as 
the  difference  between  the  integrals: 


2  (M/2)1/2  [Z1  -d(l/r1)/(x.y)1/2 


/  e  -d(l/r)/(x-y)1/2] 

1 

r,  (60 


(2y)  Z1  d  (1/r . 


-  l/r^/fx-y) 


Klein  established  the  following  definition: 


2g(y)  =  l/r2  -  1/*! 


Differentiating  with  respect  to  y  yields: 


2dg  (y)  =  (d/dy)  (l/r2)  -  (d/dy)  (l/r^ 


Changing  the  notation  yields: 


2g'dy  =  d(l/r2)  -  d(l/rx) 


Substituting  Eq  (64)  into  Eq  (61)  yields: 


P  *  2  (2y)1/2  f  [g  ’  (y)dy/(x-y)  1//2] 


^  .VV-T4-  '  -  L-  ‘  -.v7\  u  .  -AV" 


with  the  limits  of  integration  assigned  in  the  same  manner  as  for 
Eq  (52).  This  is  the  second  integral  which  is  central  to  Klein’s 
derivation . 

In  the  next  step  of  his  development,  Klein  derives  expression 
for  r  and  r  •  /  the  classical  turning  points  for  a  harmonic 
oscillator.  Klein  recognized  that  the  integrals  for  t  and  p  (Eqs 
(52)  and  (65))  were  of  the  Abelian  type.  He  multiplied  both 
sides  by  dx/(a-x)l/2  and  integrated  both  sides  from  x  =  0  to  x  =  a 
where  a  =  (E  -  V(r))^11^,  i.e.. 


a  y=  a  x  =  ct  , 

/  tdx/ta-x)1/2  =  2(2y)1'72  /  f'dy  /  dx/ [ (a-x) (x-y) ]  7 
0  y=0  x=y  (66) 

In  eq  (66) 

f ' (y)  *  f  (a)  (67) 

and 

/a  dx/ [  (a-x)  (x-y)  ]  =  /a  dx/  (-x2+y+a-ay) 1/72  (68) 

y  y 

Eq  (68)  may  be  integrated  as  follows  (Ref  6:300): 

dx/(uv)1/72  =  2/{bd)  1//2tan'1  [  (-bduv)  1//2/bv]  ,  for  bd  >  0  (69) 


36 


»  u  =  a+bx,  and  v  =  a'+dx  and  for  the  present  integral 


a  =  a  a '  =  -y 

b  *  -1  d  =  1 


bd  =  -1 

Then  the  integral  can  be  evaluated  to  be  tt  . 
Finally,  the  left  side  of  Eq  (68)  is  evaluated  as: 

/“  dx/(a-x)1/2  =  2(2y)1/2  f(a) 


olving  Eq  (71)  for  f(a): 

f(a)  =  1/(2  tt(2  p) 1/2)  /Q  dx/(  a-x) 1/2 

0 

for  a  *  y: 

y 

f(y)  =  1/(2"(2V) 1/2)  /  dx/(y-x)1/2 

0 

Similarily,  Eq  (65)  for  g  may  be  transformed  into  the 
.owing : 

g(y)  *  1/(2  tt(2  \i)1/'2)Sy  q  dx/(y-x)1/2 

0 


From  the  definitions  for  x  given  in  Eq  (44) : 


dx  *  dE 

itituting  Eq  (74')  into  Eq  (34)  yields: 

t  dx  =  dl 

illing  that  E  is  a  function  of  I  and  <  ,  and  using  Eq  (47) 
(73)  becomes: 

f(V)  =  1/(2tt  (2p  )  1/2)  /  dl/  ( V-E  ( I  ,k  )  )  1/2 

0 

From  Eqs  (34),  (35'),  and  (74*): 

3x  »  p  dE  =  (9  E/3< )  dE  =  x  (3e/9<)  (3  I/T )  =  (3  E/3<)3  I 

Eq  (74)  becomes: 

I  • 

3  (V)  *  1/ <  2  tt(  2  p)  1/2)  /  3  [E  ( I  ,<  )  /  (3k  )  dl  ]  /  (V-E  ( I  ,<  )  )  1/2 


TABLE  10 

RKR  and  IPA  Potential  for  the  a-State  of  Lead-Ovide 


GvCcnT1)  R- ( A)  R+(A)  Gv(cm_1)  R-(A)  R+(A) 


240.13 

2.075 

2574.13 

1.928 

4785.63 

1.868 

6374.63 

Y 

*nn  - 

1.827 

1.5566X10 

240.23 

716.87 

1188.55 

1655.34 

2117.24 

2574.24 

3026.34 

3473.56 
3915.90 

4353.34 
4785.85 
5213.45 
5636.23 

6054.34 
6467.82 
6876.59 


2.054 

2.011 

1.982 

1.960 

1.942 

1.925 

1.911 

1.899 

1.887 

1.877 

1.868 

1.859 

1.850 

1.842 

1.834 

1.827 

Y00  = 


2.192 

2.250 

2.293 

2.329 

2.363 

2.392 

2.421 

2.448 

2.457 

2.501 

2.527 

2.552 

2.576 

2.600 

2.624 

2.648 

0.16759 


25. 

.163 

24, 

.988 

24. 

.822 

24, 

.658 

24. 

.495 

2 

24.157 

23.979 

23.795 

23.607 

23.419 


23.060 

22.893 

22.729 

22.562 


TABLE  9 

RKR  and  IPA  Potential  for  Ground  State  of  Lead-Oxide 


RKR 


IPA 


GvCem"1 ) 

R-(A) 

R  +  (  A) 

-  —  — -  ■  -  — 

Gv ( cm" 1 ) 

R-(A) 

R+(b 

BvCcm"1 )X100 

/2 

1 

.921 

359.62 

1.868 

1.980 

359.918 

1.867 

1.979 

30.687 

1073.52 

1.832 

2.027 

1073.901 

1.831 

2.027 

30.469 

1780.34 

1.809 

2.062 

1780.716 

1.808 

2.062 

30.256 

2480.10 

1.790 

2.092 

2480.445 

1.791 

2.092 

30.048 

3172.78 

1.775 

2.118 

3173-120 

1.776 

2.119 

29.844 

3858.39 

1.762 

2.143 

3858.743 

1.763 

2.144 

29.645 

4536.93 

1.751 

2.167 

4537.299 

1.752 

2.168 

29.452 

5208.40 

1.740 

2.189 

5208.773 

1.741 

2.190 

29.264 

5872.79 

1.731 

2.21 1 

5873.161 

1.732 

2.212 

29.082 

6530.11 

1.722 

2.231 

6530.471 

1 .723 

2.236 

28.905 

7180.36 

1.714 

2.251 

7180.713 

1.715 

2.253 

28.731 

7823.54 

1.706 

2.271 

7823.891 

1.707 

2.273 

28.557 

8259.65 

1.699 

2.291 

8459.997 

1 .700 

2.292 

28.378 

9088.68 

1.693 

2.311 

9089.028 

1.693 

2.311 

28.192 

9710.64 

1.686 

2.330 

9711.015 

1.686 

2.330 

27.998 

10325.53 

1.680 

2.349 

10326.059  . 

1.680 

2.349 

27.796 

Y 

*00 

=  0.0233501 

n 

o 

o 

>* 

0.37952 

nstants 

State 

X 

a 

A 

Te* 

0.0 

16024.9(1.45) 

19862.6(1.5) 

*10* 

720.97(0.36) 

481.5(0.70) 

443.3(0.8) 

Y20* 

-  3.536(0.025) 

-  2.45(0.07) 

-  0.54(0.24) 

*01 

0.307519** 

0.252(0.10)* 

0.2588** 

Y11 

-  1.9167X10-3** 

-  1.6761X10**3# 

1.4X10"3** 

*02 

2.2X10"7** 

2.2X10_7## 

3.3X10'7** 

Mass  :  Oxygen  -  1.5994  amu;  Lead  -  207.19  amu. 
Source:  Linton  and  Broida  (Ref  33)  . 

Source:  Suchard  (Ref  45) . 


■Source:  Calculated  from  Linton  and  Broida' s  using  Eq  (104). 

Source:  Used  X-state  value  for  a-state  Yq2'  T^e  a~stat®  Yg2  be 

calculated  according  to  Eq  (104').  The  error  due  to  this 
substitution  should  be  small. 


Lead-Oxide  RKR-IPA  Curves 

To  demonstrate  the  use  of  the  RKR-IPA  program,  the  x-,  a-, 
and  A-states  were  selected  for  evaluation.  Curves  were  drawn  for 
the  x-,  a-,  and  A-states  using  the  results  of  the  RKR  routine 
alone  or  using  the  combined  RKR  and  IPA  routine  (RKR-IPA  curve) . 
The  constants  used  to  generate  these  curves  for  each  state  are 
presented  in  Table  8.  The  RKR  and  IPA  results  are  presented  for 
the  X-  and  a-states  in  Tables  9  and  10,  respectively.  The  RKR 
data  for  the  A-state  is  presented  in  Table  11.  The  resulting 
potential  energy  curves  are  presented  in  Figure  2  for  the  X-,  A-, 
and  the  a-states.  For  the  A-state,  the  curve  presented  in  Figure 
2  is  an  RKR  curve,  while  for  the  X-  and  a-states,  they  are  RKR- 
IPA  curves. 

The  ultimate  use  of  a  potential  should  be  to  predict  what 
transitions  should  occur.  A  sophisticated  analysis  would  involve 
the  calculation  of  Franck-Condon  factors  from  the  wave  functions 
resulting  from  the  solution  of  the  Schroedinger  wave  equation. 

But  it  is  possible  to  arrive  at  certain  conclusions  without 
calculation  of  the  Franck-Condon  factors.  For  example, 
examination  of  the  a-state  and  X-state  curves  show  that  no 
overlap  exists  between  the  v'=*2  energy  level  for  the  a-state,  and 
v"*0  for  the  X-state.  This  indicates  that  such  a  transition  is 
classically  forbidden.  For  example,  a  transition  observed  at 
15898  cm-^  had  been  assigned  to  the  a-X  (0,0)  transition. 


Reevaluation  showed  that  the  15898  band  head  should  be  the  a-X 


part  of  the  present  program  but  are  here  briefly  discussed.  One 
approach  Vidal  used  for  quasibound  states  involved  starting  the 
integration  of  the  SWE  at  small  internuclear  distances  and 
looking  for  the  maximum  of  the  internal  amplitude  inside  the 
centrifugal  barrier  using  a  Breit-Wigner  parametr ization ) # 
Vidal  also  used  a  second  approach  in  the  same  work(^).  For  the 
quasibound  states,  he  introduced  an  artificial  barrier  at  large 
internuclear  distances  permitting  the  use  of  the  Numerov-Cooley 
method.  The  eigenvalues  he  found  in  this  manner  were  slightly 
higher  than  the  energy  eigenvalues  derived  from  the  maximum  of 
the  internal  amplitude.  Proper  choice  of  the  barrier  kept  the 
differences  within  the  standard  errors  of  the  measurement.  Using 
this  technique,  Vidal  stated  the  same  numerical  method  can  be 
used  for  both  quasibound  and  bound  states. 

Continuing  with  the  program,  as  presented  in  this  paper,  it 
is  necessary  to  transform  the  energy  used  in  the  SWE  to  units 
which  are  consistent  with  the  units  of  length  being  used,  Bohr 
radii.  Thus  far,  energy  has  been  expressed  in  terms  of  inverse 
centimeters,  wavenumbers.  To  convert  wavenumbers  to  energy  units 
consistent  with  Bohr  radii,  it  is  necessary  to  divide  by  this 
factor  . 


*Na  60.19972628 
*—  =  ■  wavenumbers 

47rca0UA  UA 


(107) 


where  a 


0.52197706  and  NA  =  Avogadro's  number. 


the  outer  turning  points  of  the  rotationless  potential.  To  avoid 
this,  Vidal  chose  a  nonlinear  interpolation  given  by: 


(r“re) ^rmax“rmin^ 


(rmax+rmin) ^re+r^ "^rraaxrmin"^rer 


(105) 


This  relation  assumes  X=1,  for  r=r  and  X=0,  for  r=rmax  and  X=1 


for  r=rmin  * 


The  interpolation  becomes  linear  for: 


re  = 


rmax+rmin 


(106) 


This  formulation  treats  the  inner  and  outer  turning  points  with 
comparable  weight  and  reduces  the  number  of  Legendre  polynomials 
in  Eq  (105). 

As  noted,  the  first  use  of  the  perturbation  technique  to 
adjust  approximate  potential  energy  curves  is  attributed  to  Hinze 
and  Kosman^6)#  The  technique  used  in  this  program  to  solve  the 
Schroedinger  equation  is  the  Nuraerov-Cooley  method  » *^) .  The 
Numerov-Cooley  method  of  solving  the  SWE  was  used  by  J.  K. 

Cashion  to  test  the  validity  of  approximate  eigenvalue  equations 
developed  by  Perkeris  for  a  rotating  Morse  oscillator ^ . 

The  routine  presented  in  this  paper  solves  the  SWE  for  bound 
states.  Vidal  discusses  techniques  for  handling  quasibound 
states.  The  techniques  for  evaluating  quasibound  states  are  not 


’.*  *  *  •  v *'»  .  • 


to  extend  calculations  beyond  RMII  and  RMAA,  as  specified  by  Eyj, 
produced  large  oscillations.  He  finally  settled  upon  a 
combination  of  Legendre  polynomials  and  an  exponential  function 
to  dampen  the  oscillations.  This  combination  improved  the 
convergence  of  the  IPA  method.  The  expression  Vidal  used  is  as 

follows(  . 


VQ(r)  =  E  ciPi(x)exp(-x2n)  (103) 

where  the  typical  range  for  n  is  1<n<5.  The  Gaussian  part  of  Eq 
(103)  provides  a  smooth  cutoff  avoiding  unphysical  oscillations. 
The  Legendre  polynomials  are  calculated  using  the  standard 
recursion  relations ^°) . 

(n+1>pn+i(x)  =  (2n+1)XPn(x)-nPn_1(x)  (104) 

where  P-|  =  1  and  ?2  ~ 

Kosman  and  Hinze^S)  used  a  linear  relationship  between  r  and 
X  such  that  X=1  for  r=4max  and  X=-1  for  r-rrain*  Vidal  found  this 
relationship  to  provide  poor  convergence  when  dealing  with 
vibrational  levels  all  the  way  to  the  dissociation  limit  and  when 
dealing  with  anharraonic  potentials.  He  stated  that  the  reason 
for  the  poor  convergence  in  the  case  of  a  highly  anharmonic 
potential  is  that  a  linear  interpolation  tends  to  optimize  only 
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v0(r),  being  sought  may  be  expressed  as  a  sum  of  the  approximate 
potential  and  some  delta  potential. 

V0(r)  =  V0(r)  +  AV0(r)  (99) 

In  turn,  it  follows  that  the  true  Hamiltonian  differs  from  the 
approximate  zero  order  Hamiltonian  (HQ  +  Hrofc)  by  AV  (r) : 

H0  +  Hrot  "  H0  +  Hrot  +  Av0(r)  (100) 

First  order  perturbation  theory  gives  the  following  relationship 

between  the  perturbation  to  the  Hamiltonian  AV0  (r)  and  the 
energy  change  AEvjJ 

AEvJ  =  <^j  AV0(r)  4j>  (10D 

To  obtain  a  correction  to  V(r),  a  mathematical  expression  of 
the  following  form  is  assumed  to  represent  AVg  (r) : 

Av0  (t)  =  jc^U)  (102) 

The  selection  of  the  specific  form  of  the  function  f^(r)  is 
critical  to  the  convergence  and  solution  of  the  SWE. 

Kosman  and  Hinze  chose  Legendre  polynomials  (x)  to 
represent  the  f^(t)  f unctions ( 36) .  Vidal  found  that  attempting 


This  formulation  has  been  simplified  by  neglecting  terms 
describing  electronic  coupling  contributions  . 

The  rotational  motion  of  the  molecule  is  described  by: 


(97) 


In  the  inverted  perturbation  approach,  one  starts  with  an 
approximate  potential  Vq(>*),  for  example  an  RKR  potential.  Eq 
(95)  is  solved  numerically  for  the  zeroth  order  eigenvalues  EyJ 

and  the  radial  wavefunetion  <|>vj(r)* 

Then  an  energy  correction  AEyJ  is  calculated  according  to  the 
following  formula: 

iEvJ  =  EVJ  -  EW  (98) 

where  Evj  is  the  measured  terra  value  as  calculated  from  the 
spectroscopic  constants  and  quantum  numbers  v  and  J.  If  the 
difference  between  the  experimentally  determined  term  values  and 
the  eigenvalues  obtained  from  the  SWE  is  sufficiently  small,  the 
calculation  is  stopped. 

If  the  difference  exceeds  a  specified  limit,  then  one 


proceeds.  As  stated,  the  goal  is  to  obtain  some  correction  to 
V(r),  i.e.,  AV(r),  such  that  the  calculated  and  experimental 
EvI's  agree.  Following  this  reasoning,  then  the  true  potential, 


with  the  SWE.  This  leads  to  the  use  of  the  inverted  perturbation 
approach.  The  goal  of  the  inverted  perturbation  approach  (IPA) 
is  to  adjust  the  potential  energy  curve  of  rotationless 
molecules,  Vg(*)/  so  that  the  quantum  mechanical  eigenvalues, 

Evj,  obtained  from  a  solution  of  the  Schroedinger  wave  equation 
agree  in  a  least-squares  sense  with  the  measured  term  values, 
T(v,J) . 

The  IPA  technique  was  first  demonstrated  by  Hinze  and 
Kosman^36).  Vidal  and  Scheingraber expanded  the  use  of  the 
IPA  method.  Vidal ( 37 ' 38 ) ,  in  conjunction  with  several  authors, 
has  applied  the  IPA  method  to  a  number  of  molecules. 

The  Schroedinger  wave  equation  (SWE)  for  a  vibrating  rotator 
(Eq  (4))  is  the  basis  of  the  IPA  method.  To  develop  the  IPA 
method,  the  SWE  can  be  expressed  as  follows 


(Ho  +  Hrot)^vj(r)  «  EvJ*vJ<r>  (95) 

E°vJ  is  the  energy  eigenvalue  specified  by  the  vibration  and 
rotation  quantum  numbers.  H0  is  the  Hamiltonian  of  the 
nonrotating  molecules  and  is  made  up  of  the  terms.: 


H 


0 


dr' 


+  V0(r) 


(96) 


where  V0(r)  is  the  potential  energy  of  the  rotationless  molecule. 
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Using  these  definitions,  Eq  (86)  for  g  becomes:'1-3  ‘ 


g ( v )  =  /  B(v')(G(v)-G(v'))"1/2dv»  (93) 

v0 

From  the  definitions  in  Eqs  (22)  and  (62)  for  f  and  g  the 
following  expression  for  inner  and  outer  turning  points  may  be 
obtained : 


r+  =  (f2  +  f/g)1/2  +  f 


(9*0 


Solutions  of  the  equations  for  f  and  g  are  discussed  in  the 
section  on  the  RKR  Program  by  C.  R.  Vidal'7.  The  Dunham 
coefficients  calculated  from  spectroscopic  constants  (i.e.,  from 
Program  DUNCON)  are  input  to  the  program,  program  RIPA.  By  use 
of  this  program  the  integrations  are  performed  by  means  of  a 
Gaussian  integration  routine  and  an  RKR  potential  energy  curve  is 
produced  by  means  of  a  least  squares  fit  to  the  inner  and  outer 
portions  of  the  curve.  These  calculations  are  performed 
separately  for  each  electronic  state  of  interest. 

Inverted  Perturbation  Analysis 

Having  constructed  a  seraiclassical  RKR  curve,  the  next  step 
is  to  determine  if  that  curve  is  consistent  with  the  Schroedinger 
wave  equation;  or,  even  better,  to  adjust  the  curve  so  it  agrees 


Eq  (85)  becomes: 


v ' 

f  *  /  (G  (v)  -  G  (v ' )  )  *"1/2  dv’  (88) 

v0 

where  Kaiser  defined  the  lower  limit  of  integration  by  the 

following : 


G(v0) 


"  -Y00 


(89) 


where  Y00  is  the  Dunham  constant  previously  defined  in  Eq  (6) . 
The  variable  v0  j.n  G(v0)  may  be  expressed  in  terms  of  Dunham 
coefficients  as: 


v0  =  -1/2  -  A  =  -1/2  (Y00/Y10)(l  +  Y00Y20/Y10  +  •••)  (90) 


For  J 


0 


(34) 


3E/3J(J+1)  =  nZ0  Yn.(v+l/2)n  =  B (v )  (91) 

wnere  B(v)  is  the  spectroscopic  term  defined  as  follows: (13^ 

Bv  *  Be  “  ot  (v+1/2)  +  y(v+l/2)2... 


(92) 


3/  3<  =  2jj/h‘t(3/3  J(J+1)) 


Finally,  the  expressions  for  f  and  g  become: 


f  (V,  J  (J+l) )  *  h/(2y)1/2 


v ’ +1/2  d  (v  +  1/2) 


(V-E (v+1/2, J (J+l) ) ] 


g  (V,  J  (J+l)  =  (2y)1/2/h 


v '  +1/2  3E/3(J(J+1)  (v+1/2) 


d (v+1/2) 


(V-E (v+1/2, J (J+l) ) ] 


where  the  upper  limit  of  integration  v'  is  selected  such  that 
V  =  E (v ' +1/2 , J ( J+l) )  for  the  fixed  value  of  J(J+1)  used  in  the 
integrals  . 

If  the  RKR  calculations  are  evaluated  for  J=0  then  Eq  (2) 
reduces  to: 


T  "  i=l  Yi0  (v+1/2)1  =  G(v) 


more  commonly  referred  to  as 


As  E  ( I ,  <)  approaches  V  in  Eqs  (76)  and  (78),  the  integrals 
become  infinite.  To  avoid  these  singularities,  Klein  instead 
evaluated  the  following  expression: 

S(V,<)  =  l/(ir(2y) I/2) /(V-E(I,  <)  )  1/2dl  (79) 

This  expression  has  the  following  relation  to  the  f  and  g  of 
Eqs  (76)  and  (78): 

f  =  3S/3V  ;  g  ■  -3S/3<  (80) 

As  shown  in  Eq  (2),  spectroscopic  energy  levels  can  be 
represented  as  power  series  of  (v+1/2)  and  J(J+1).  It  is 
convenient  then  to  express  the  integrals  for  f  and  g  in  terms  of 
the  quantum  numbers  v  and  J. 

The  expression  for  the  radial  action  variable  in  quantum 
mechanical  terms  is 

I  *  (v+1/2) h  (81) 

Then 

d ( I )  *  d  (v+1/2) h  (82) 

The  quantum  mechanical  expression  for  k  is: 


<  *  P2/2  «  J(J+1)  h2/2 


(83) 


TABLE  11 


(3,2)  transition.  Both  assignments  were  acceptable  within  the 
error  limits  of  fit.  The  relative  positioning  of  the  potentials 
indicate  that  only  the  (3,2)  transition  is  possible.  The 
calculated  value  for  the  a-X  (3,2)  transition  is  15898.4  cm"^-. 

The  next  step  should  be  the  generation  of  a  new  set  of 
constants  using  the  formulas  for  Gv  and  Bv  and  the  newly  obtained 
energy  eigenvalues  Gv  and  Bv  values  generated  by  the  IPA  routine. 
This  technique  was  applied  to  the  ground  state  IPA  potential  of 
Table  9.  From  this,  the  following  constants  were  generated  for 
the  ground  state: 

*10  =  720.991(0.058)  cm'1 

Y20  =  3.5357(0.00479) 

The  calculation  of  new  values  for  the  rotational  constants  Y0^f 
Yn,  etc.,  from  Bv  was  not  performed. 

The  calculations  using  only  Vidal's  RKR-IPA  method  were 
halted  at  this  point.  Attempts  to  extend  the  curves  to  the 
dissociation  limit  in  one  jump  produced  a  potential  which 
diverged  as  the  curve  approached  the  dissociation  limit  (see 
Figure  3) .  Not  having  data  for  vibrational  levels  above  about 
v=15,  the  iterative  techniques  of  gradually  extending  the  curve, 
as  discussed  by  Vidal  and  Stawalley  (Ref  43),  were  not  attempted. 
They  used  the  IPA  method  to  improve  the  RKR  curve  at  lower 
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-  rkr  Curve  (Linton  and  Broida,  Suchard  Constants) 

- -  "Extended  Constants,”  RKR- I PA  Curve 

Figure  3.  Lead-Oxide  X-State  Potential  Energy  Curves 


levels,  and  then  on  successive  trials,  extended  the  curve  slowly 
to  the  dissociation  limit. 

The  reader  is  referred  to  Vidal's  works  listed  in  the 
Bibliography  for  further  information  concerning  the  generation  of 
new  constants  and  the  extension  of  potential  curves  to  the 
dissociation  limit  using  the  IPA  routine. 

Methods  for  Extending  RKR  Curves 

In  many  cases,  because  of  the  lack  of  experimental  data,  the 
constants  necessary  to  construct  an  RKR  curve  may  not  be 
obtainable  from  straightforward  least-square  fits  or  their 
quality  may  be  insufficient  to  construct  accurate  curves.  Two 
techniques  for  improving  existing  constants  which  may,  in  turn, 
be  used  to  produced  improved  RKR  curves  are  presented  here^22). 

The  potential  energy  curve  of  a  diatomic  molecule  should  have 
a  minimum  at  some  intermediate  distance,  approach  the 
dissociation  limit  for  large  internuclear  separations  and 
increase  rapidly  as  the  atoms  approach  each  other.  Ideally, 
curves  generated  using  the  RKR  method  should  satisfy  these 
criteria. 

For  many  diatomic  molecules,  only  two  and,  at  most,  three 
vibrational  constants  are  available.  The  available  rotational 
constants  are  even  more  limited.  Attempts  to  generate  curves 
based  upon  these  constants  using  the  RKR  method  can  be  expected 
to  produce  curves  which  do  not  meet  the  dissociation  criteria  or 
which  do  not  exhibit  an  increasing  potential  for  small 
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internuclear  distances.  These  results  occur  because  the 


experimental  data  upon  which  the  constants  are  based  may 
represent  only  a  fraction  of  the  vibrational  levels  which  exist 
between  the  bottom  of  the  potential  well  and  the  dissociation 
limit. 

To  develop  potential  energy  curves  which  obey  the  known 
constraints,  Leroy  and  Burns  have  reported  on  a  new  technique 
which  shows  promising  results.  While  they  start  with  constants 
obtained  from  experimentally  observed  transitions,  they  then 
adjust  these  constants  slightly  until  the  new  RKR  potential 
energy  curve  is  most  consistent  with  all  of  the  known 
constraints. 

Leroy  and  Burns  used  their  knowledge  of  the  desired  shape  to 
adjust  the  constants.  They  worked  with  the  Gv  and  Bv  quantities 
as  expressed  in  Eqs  (87)  and  (91) .  Expressed  in  terms  of  these 
two  quantities,  their  criteria  for  adjusting  the  molecular 
constants  are: 

1.  V(  )  assymptotically  approaches  the  dissociation  energy 
(De)  and  the  difference  between  energy  levels,  Gv+ 1/2,  becomes  0 
for  the  same  value  of  v.  Gv+l/2  is  defined  as: 

Gv+1/2  "  G(v  +  D  -  G<v>  <108) 


2.  The  slope  of  the  inner  portion  of  the  RKR  curve  must  be 
negative.  Hence: 


dV(r)/dr  <  0 


(109) 


3.  The  slope  of  the  inner  portion  of  the  curve  must  become 
increasingly  steeper;  that  is,  the  second  derivative  with  respect 
to  the  internuclear  separation  must  be  positive: 

d2V(r)/dr2  >  0  (110) 

A  study  of  Eq  (88)  reveals  that  as  gw+i/2  becomes  smaller, 
the  quantity  f  becomes  larger.  In  turn,  Eq  (93)  shows  that  f 
determines  the  width  of  the  RKR  curve.  Examination  of  Eq  (93) 
further  reveals  that  the  value  g  determines  the  center  of  the 
potential  curve  for  a  given  energy  level.  In  turn,  Eq  (92)  shows 
that  having  determined  the  constants  necessary  to  generate  the 
G(v)'s,  the  other  factor  affecting  the  value  of  g  is  B(v),  a 
quantity  whose  value  is  determined  by  the  rotational  constants, 
*ni.  From  Eqs  (88),  (92),  and  (93),  it  follows  that  the  width  of 
the  curve,  as  determined  by  the  RKR  method,  and  the  value  at 
which  gv+i/2  =  0  depends  only  upon  the  vibrational  constants 
which  make  up  Gv.  On  the  other  hand,  having  selected  suitable 
constants  for  Gv,  the  behavior  of  g,  the  rate  at  which  the  outer 
portion  of  the  curve  approaches  D er  and  the  behavior  of  the  slope 
of  the  inner  portion  of  the  curve  depends  solely  upon  the 
constants  which  constitute  By- 

Following  this  logic,  Leroy  applied  the  criteria  in  1  above 
first.  If  the  constants  obtained  from  experimental  data  do  not 
satisfy  this  criteria,  the  value  of  the  last  experimental 


constant  is  adjusted  or  a  value  is  selected  for  the  next 
constant,  Yn0/  in  the  series  given  by  expression  (8) .  The 
contribution  from  that  constant,  Yn0(v+l/2)n»  is  then  subtracted 
from  the  G(v)  values  calculated  from  the  original  experimentally 
based  constants.  A  least-squares  fit  is  then  performed  on  the 
adjusted  experimental  data,  Gv j : 

n-1 

GVj  =  i-  Yi0  (V3  +  1/2)1  "  *n0  (Vj  +  1/2)0  (111) 

For  a  given  YR0  the  experimental  data  fixes  the  values  of  the 
other  constants.  This  new  set  of  constants  is  tested  to 
determine  if  Gv  approaches  a  maximum  value  of  D0.  Leroy  and 
Burns  say  that  if  this  value  approaches  De  within  <2  cm-1  then 
the  criteria  in  1  is  satisfied.  If  1  is  not  satisfied,  then  a 
new  trial  value  for  Yn0  is  chosen  and  the  process  repeated  until 
1  is  satisfied. 

At  this  point,  a  new  RKR  curve  is  generated  using  the  new 
vibrational  constants,  Y^0*  and  the  experimental  rotational 
constants,  Y^.  The  curve  is  then  evaluated  according  to 
criteria  (2)  and  (3).  If  the  inner  portion  of  the  curve 
diverges,  the  process  followed  for  adjusting  the  constants  for  Gv 
is  repeated  for  the  constants  making  up  Bv.  The  value  for  the 
last  experimental  constant  is  adjusted  or  a  value  for  the  next 
rotational  constant  in  the  Yji  series  is  selected.  Then  B., 


values  given  by  the  experimental  constants  are  adjusted  as 
follows: 

Bvi  =  V  Yil(vi+l/2)i  -  [Ynl(Vi+l/2)n]  (112) 

i=0 

Then  a  least-squares  fit  is  performed  to  obtain  a  new  set  of 
rotational  constants  Yg^  through  Ynl.  The  new  constants  are  then 
used  to  produce  a  new  RKR  curve.  If  the  curve  does  not  satisfy 
the  criteria  in  2  and  3/  then  the  Yn^  constant  is  adjusted  and 
the  process  repeated  until  2  and  3  are  satisfied. 

In  this  manner,  Leroy  and  Burns  obtained  a  set  of  constants 
consistent  with  the  experimental  data  and  with  criteria  1-3. 

Tellinghuisen  and  Henderson  describe  a  technique  for 
constructing  RKR  based  curves  when  experimental  data  does  not 
permit  the  direct  calculation  of  the  rotation  constants . 

Their  technique  is  based  upon  a  combination  of  the  Morse  and  RKR 
potentials;  hence,  Morse-RKR  curves. 

As  explained,  vibrational  constants  provide  all  the 
information  necessary  to  calculate  f,  Eq  (88)  and,  thus,  the 
width  of  the  potential  curve,  2f,  as  defined  by  Klein  for  a  given 
vibrational  energy  level.  Having  established  the  width  of  a 
potential,  the  remaining  piece  of  information  required  to 
construct  a  potential  curve  is  the  location  of  either  the  inner 
or  outer  branches  of  the  potential  curve,  or  the  center  of  the 
potential  well,  i.e.,  g  as  obtained  from  the  RKR  calculations. 


Tellinghuisen  and  Henderson  point  out  that  inner  turning 
points  on  a  potential  curve  might  be  approximated  by  the  Morse 
potential  m 


U(r)  =  De  (i  _  e'8(r_re  })2 


where: 


D  =  dissociation  energy,  cm" 


re  =  equilibrium  internuclear  distance. 


*  *  0.12777  we(VDe)1/2 


where 


oje  »  the  vibration  frequency,  cm" 


y  a  reduced  mass,  amu 


Or  because  of  the  relation 


D  «  we2/4(wexe) 


S  may  be  expressed  as  follows: 


B  *  0.243555  (yu  x  ) 1/2 


J*-  A.  J-  J  -*  -*  — ■  *  V*  'J  •  Ji  •  V  *_•  «*■•*-  -  •  •  .  '  .  "  .  >  ! 
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Using  these  relations,  Tellinghuisen  investigated  25 
different  well-known  potentials  and  found  errors  at  the 


dissociation  limit  for  the  inner  portion  of  the  curve  to  be 
typically  less  than  0.02  a  and  re  t0  he  i-n  error  hy  seldom  more 
than  1%. 

Further,  Tellinghuisen  states  that  a  potential  curve  in  error 
bv  Ar  •  has  identical  classical  and  almost  identical  quantum 
eigenvalues  for  J  =  0.  The  wavef unctions  derived  from  a  shifted 
curve  would  be  skewed  with  respect  to  a  proper  wavefunction. 

If  possible,  Tellinghuisen  recommends  the  construction  of 
curves  by  formulas  (115)  and  (116)  using  the  experimental  and 
wexe.  If  only  e  and  De  are  known,  then  a  value  of  De  40%  larger 
than  the  experimental  De  should  be  used. 

Tellinghuisen  also  speculates  that  ae  might  be  calculated 
using  an  expression  derived  by  Pekeris  from  the  Morse 
function: 


Tellinghuisen  found  that  for  23  of  the  25  cases  he  investigated, 
the  ae  calculated  per  formula  (117)  had  an  average  absolute  error 
of  13%  and  an  average  signed  error  of  6%. 


If  one  is  seeking  the  rotational  constant  De  (Yq2),  the 

(13) 

following  relation  might  be  used: 


Tellinghuisen  cautions  against  using  Bv  values  made  up  of 
only  two  terms,  i.e., 


Bv  =  Be  _a  (v+1/2)  (119) 

because  the  Morse  expression  for  Bv  does  not  terminate  with  two 
terms  as  does  Gv.  Further,  for  high  values  of  v,  the  Morse  Bv 
will  differ  from  the  Morse-RKR  Bv.  The  Morse-RKR  Bv  value  must 
be  determined  by  numerical  methods. 

b)  Extension  of  the  X-State  Lead-Oxide  Potential  Energy  Curve  to 

the  Dissociation  Limit, 

_ ’  e 

The  potential  energy  curve  for  the  ground  state  of  lead-oxide 
wa3  selected  as  the  likely  candidate  for  extension  to  the 
dissociation  energy,  Dg.  This  selection  was  based  upon  the  fact 
that  the  constants  for  the  ground  state  have  the  smallest 
standard  errors. 

The  vibrational  constants  of  Linton  and  Broida  were  selected 


as  starting  values  because  of  their  small  error  limits.  The 
creation  of  potential  energy  curves  also  requires  the  knowledge 


tational  constants.  The  rotational  constants  were  selected 
Suchard^4^.  Table  12  contains  the  original  constants  and 
inal  set  of  constants  used  in  extending  the  curve  to  the 
iciation  limit. 

'he  method  prescribed  by  Leroy  was  used.  Linton  and  Broida's 
:ants  were  used  to  generate  a  set  of  G(v)'s  according  to  Eq 
for  v=0  to  15.  A  trial  value  was  selected  for  Y^0  and  the 
c  constants,  Y30,  Y20,  and  Y10  were  obtained  by  performing  a 
t-squares  fit  to  the  value  obtained  from  Eq  (111).  Program 
ON  was  modified  to  perform  the  fit.  The  process  was  repeated 
1  the  criterion  of  <»v+1/2  *  0  at  Dfi  was  satisfied. 

Suchard's  constants  were  used  to  generate  a  set  of  Bv's  for 
to  15  according  to  Eq  (92)  .  A  trial  Y^  was  selected  and  new 
and  Yj^  constants  were  obtained  from  a  least  squares  fit  to 
values  calculated  by  Eq  (112)  . 

The  new  vibrational  constants  and  rotational  constants  were 
as  the  inputs  for  the  RKR-IPA  program.  The  process  of 
sting  the  rotational  constants  was  repeated  until  the 
eria  expressed  in  Eqs  (108),  (109),  and  (110)  were  satisfied. 
The  last  set  of  constants  calculated  are  presented  in  Table 
s  the  extended  constants.  The  Final  Set  of  Turning  Points, 
utput  by  the  IPA  routine  was  used  to  create  the  potential 
e  (solid  line)  in  Fig  3  which  satisfies  the  dissociation 
eria.  The  dashed  line  in  Fig  3  is  the  RKR  curve,  as  defined 
he  original  constants  of  Linton  and  Broida,  and  Suchard.  It 


TABLE  12 


Constants  Used  to  Extend  the  Lead-Oxide  X-State 
Potential  Energy  Curve  to  the  Dissociation  Limit,  Dq 


Constants 

Original  Values 

Extended  Values 

Y10 

720.97(0.36) 

721.062849 

y20 

-  3.536 

-  3.57101904 

Y30 

3.85076374XI0"3 

y40 

-  1.29X10"4 

Y01 

0.307519 

0.30726775 

Y11 

-  1.9167X10-3 

-  1 . 7667X1 O”3 

Y21 

-  1.5X10"5 

y02 

0.22X10-6 

0 . 22X1 0"6 

De  =  31570+410.  era"1 

should  be  noted  that  the  constants,  as  presented  in  Table  12, 
will  only  produce  a  curve  in  Fig  3  when  input  into  an  RKR-IPA 
routine.  Examination  of  the  Final  Set  of  Turning  Points  yields  a 
minimum  inner  turning  radius  of  about  1.62  angstroms  for 
vibrational  quantum  numbers  from  v=33  to  v=69.  Examination  of 
the  turning  points  for  the  initial  potential  generated  by  RKR, 
listed  at  the  beginning  of  the  RKR-IPA  program  output  shows  RMIN 
turning  points  as  small  as  1.54  angstroms  for  v=74.  For  v=69, 
the  last  vibrational  quantum  number  investigated  by  the  IPA 
routine,  the  RKR  program  returned  an  RMIN  of  1.59  angstroms. 

This  decreasing  value  of  RMIN  indicates  that  the  left  branch  of 
the  curve  would  diverge  just  as  shown  by  the  dashed  line  on  the 
left-hand  side  in  Fig  3.  The  turning  points  of  the  RKR  potential 
were  of  sufficient  quality  that  the  IPA  program  was  able  to 
correct  the  potential  to  yield  the  solid  curve  of  Fig  3.  In  this 
work,  the  procedure  described  above,  was  the  final  one  used; 
however,  continued  refinements  are  possible  by  the  use  of  the 
RKR-IPA  program.  The  next  step  should  be  to  take  the 
G(v)  [U ( R ) ]  and  B(v)  [BV*100]  values  as  output  by  the  IPA  program 
and  perform  a  least-squares  fit  to  obtain  a  new  set  of  YnQ  and 
Ynl  constants.  These,  in  turn,  would  be  used  as  input  for  the 
constants  from  the  experimental  data  and  to  calculate  the  higher 

Yn+,  v+  constants. 

The  convergence  of  the  IPA  routine  became  very  sensitive  as 
the  sought  for  value  of  the  last  vibrational  constant  was 


approached.  For  example,  the  IPA  routine  would  not  converge  when 
values  of  Y21  =  -1.25x10"5  or  y21  =  -1.75x10-5  were  used.  The 
program  would  converge  for  values  of  Y21  larger  and  smaller  than 
these,  but  the  shape  of  the  resulting  curve  diverged  from  the 
desired  shape.  An  indication  of  the  convergence  of  the  IPA 
routine  may  be  obtained  by  looking  at  the  next  to  the  last  output 
of  the  IPA-RKR  program,  "The  Summary  of  Errors  of  the  Inverted 
Perturbation  Approach.” 

Due  to  the  fact  that  the  final  set  of  turning  points 
fluctuates  about  1.62  angstroms  and  the  potential  is  in  reasing 
asymptotically  along  1.62  angstroms,  multiple  potential  energies 
are  presented  for  the  same  RMIN  turning  point  value.  Thus,  when 
searches  are  performed  using  a  turning  point  radius  value  [XXX] 
as  in  DO-loop  55  of  POTTAB,  and  PLYNN  called  from  this  loop,  the 
routine  stops  at  the  first  value  satisfying  the  requirement  [XXX] 
X(J)  (Subroutine  PLYNN).  The  end  result  is  that  potential 
energies  corresponding  to  turning  points  of  magnitude  of  less 
than  about  1.62  angstroms  in  the  final  potential  should  not  be 
accepted  as  valid  without  close  inspection.  Improved  constants 
and  successive  iterations  of  the  RKR-IPA  routine  may  remove  this 
discrepancy . 
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